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ABSTRACT 

We study the properties, evolution and formation mechanisms of isolated field ellip- 
tical galaxies. We create a 'mock' catalogue of isolated field elliptical galaxies from 
the Millennium Simulation Galaxy Catalogue, and trace their merging histories. The 
formation, identity and assembly redshifts of simulated isolated and non-isolated el- 
liptical galaxies are studied and compared. Observational and numerical data are used 
to compare age, mass, and the colour-magnitude relation. Our results, based on sim- 
ulation data, show that almost seven per cent of all elliptical galaxies brighter than 
— 19 mag in B— band can be classified as isolated field elliptical galaxies. Results also 
show that isolated elliptical galaxies have a rather flat luminosity function; a number 
density of - 3 x 10 _6 /i 3 Mpc~ 3 mag , throughout their B— band magnitudes. Iso- 
lated field elliptical galaxies show bluer colours than non-isolated elliptical galaxies 
and they appear younger, in a statistical sense, according to their mass weighted age. 
Isolated field elliptical galaxies also form and assemble at lower redshifts compared to 
non-isolated elliptical galaxies. About 46 per cent of isolated field elliptical galaxies 
have undergone at least one major merging event in their formation history, while the 
same fraction is only ~ 33 per cent for non-isolated ellipticals. Almost all (~ 98 per 
cent) isolated elliptical galaxies show merging activity during their evolution, pointing 
towards the importance of mergers in the formation of isolated field elliptical galaxies. 
The mean time of the last major merging is at z ~ 0.6 or 6 Gyrs ago for isolated 
ellipticals, while non-isolated ellipticals experience their last major merging signifi- 
cantly earlier at z ~ 1.1 or 8 Gyrs ago. After inspecting merger trees of simulated 
isolated field elliptical galaxies, we conclude that three different, yet typical formation 
mechanisms can be identified: solitude, coupling and cannibalism. Our results also 
predict a previously unobserved population of blue, dim and light galaxies that fulfill 
observational criteria to be classified as isolated field elliptical galaxies. This separate 
population comprises ~ 26 per cent of all IfEs. 

Key words: methods: numerical - galaxies: formation - galaxies: evolution - galaxies: 
ellipticals - galaxies: isolated - large-scale structure of Universe. 



1 INTRODUCTION 

Galaxies reside in three different environments: clusters, 
groups and voids. The majority, more than 50 per cent, of 
galax ie s are found in gr o ups an d clusters Jrlumason et al] 
19561: iHuchra fe Gelid Il982l; iGeller fc Huchral Il983t 



Nolthenius fc White! Il987l ; iRamella et alj l2002h . This 



particularly true for elliptical galaxies, as t he morphology - 
density relation has shown (|Oemlerl Il974l : [Dressier! !l98C0. 



The merger hypothesis (|Toomre fc Toomrd Il97j ) suggests 
that the product of the merger of two spiral galax- 
ies will be an elliptical galaxy. If this holds, then the 
probability to find an elliptical galaxy grows in environ- 
ments with high densities and low velocity dispersions, 
e.g. groups of galaxies. Therefore, local environment is 
thought to play a crucial role in g alaxy f ormation and 
evolu t ion (e.g. iFarouki fc Shapirol Il98ll: iBalogh et alj 



19991: Kauffmann et alj 19991 : Lemson fc Kauffmann 19991 : 
Moore et alj|l999h . 
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In observations it is relatively easy to study cluster and 
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group galaxies. Because of this and the potential formation 
mechanism of elliptical galaxies, it is not surprising that 
most studied ellipticals are found in clusters and groups. 
Therefore, the detailed properties of isolated field elliptical 
(IfE) galaxies have not been extensively studied or very well 
understood; there are only a few observational studies and 
the surveys are small. Moreover, the formation mechanisms 
and evolutionary paths of these "lonely" elliptical galaxies 
are not yet well understood. 

In the past two decades several observational projects 
have identified and st udied the properti e s of isolated field 
elliptical galaxies (e.g. | Smith et al] 12004: iReda et alJ 12004 : 
iHernandez-Toledo et alj|2008l : iNorberg et alj|2008l . and ref- 
erences therein). In many studies, based on observational 
data, different p ossible formation scenarios ha v e been pro- 
posed (see e.g. iMulchaev fc Zabludofj 1 19991 ; | Reda et al] 
l2004l2007h . ranging from clumpy collapse at an early epoch 
to multiple merging events. Also equal-mass mergers of 
two massive galaxies or collapsed groups have been sug- 
gested. Theoretical studies predict that isolated ellipticals 
are formed in relatively recent mergers of spiral galaxy 
pairs, while large isolated ellipticals ma y be the result of 
merging of a small gr oup of galaxies (e.g. IJones et al. 2000; 
iD'Onghia et~ai1l2005r i. 

Observational studies have shown that several IfEs re- 
veal a number of features such as tidal tails, dust, shells, 
discy and boxy isophotes and rapidly rota t ing discs (e.g. 
Reduzzi et all Il996l; IReda et al l I200I 120051 ; lHau fc Forbesl 
20061 : Hernandez- Toledo et al] 20081) indicating recent 
merge r and/or accretion events~ lrlernandez- Toledo et al] 
(2008) concluded that at least 78 per cent of their isolated el- 
liptical galaxies show some kind of mophological distortion, 
and suggested that thes e galaxies suffered late dry mergers. 
IReda et al] l|2004l . 120051 ) compiled a sample of 36 candidates 
of isolated early-type galaxies and studied their properties, 
and concluded that a major merger of two massive galaxies 
could explain most observed features. They also concluded 
that a collapsed poor group of a few galaxie s is a p ossible for- 
mation scenario. However, iMarcum et al] (|2004l ) studied a 
similar sample of isolated early-type galaxies, and concluded 
that isolated systems are underluminous by at least a mag- 
nitude co mpared with obje cts identified as merged group 
remnants. IReda et al] (|2007l ) concluded that mergers at dif- 
ferent redshifts of progenitors of different mass ratios and 
gas fractions are needed to reproduce the observed proper- 
ties of IfEs. 

However, some IfEs ha ve not shown any evidence of re- 
cent merging acti v ity (e .g. lAars et "all l200ll : iDenicolo et al] 
|2005| ).Ia ars et all (|200ll ). who studied a sample of nine iso- 
lated elliptical galaxies, identified five galaxies that were lo- 
cated in environments similar to those of loose groups, while 
the environments of the remaining four galaxies were con- 
firmed to be on low density. All galaxies showed smooth, 
azimuthally symmetric profiles, with no obvious indications 
of dust lanes, nascent spiral structure or star-forming re- 
gions. It is possible that the merging events have happened 
in distant past, and all s igns of merging events have been 
wiped out. iMihosI (|l995h found out, by using a combina- 
tion of numerical simulation and synthesized Hubble Space 
Telescope (HST) Wide-Field Planetary Camera 2 (WFPC2) 
images, that merger remnants appear morphologically indis- 
tinguishable from a "typical" elliptical ^ 1 Gyrs after the 



galaxies merged, while ICombes et al] (|l995T ) estimated from 
numerical simulations that the time might be even less than 
0.5 Gyrs. Indeed, these times are very short in time scales 
of galaxy evolution, making it difficult to find observational 
evidence to back up formation via merging events. On the 
other hand, it is equally possible that some or even all of 
these galaxies have initially formed in underdense regions 
and developed quietly without any major mergers. 

Observational studies of isolated field elliptical galax- 
ies use similar isolation criteria as optical studies of fossil 
groups. Both classes show a large magnitude gap between 
the first- and the second-ranked galaxy. However, there is no 
criterion for IfEs that would require a presence of extended 
X-ray emission, as in case of fossil groups. Additionally, fossil 
groups are not necessarily found in low-density environments 
l|von Benda-Beckmann et alj|200sl ) like IfEs. Despite the dif- 
ferences, it is possible that the formation mechanisms of the 
two systems are similar or related to each other. Therefore, 
it is interesting to compare isolated field elliptical galaxies 
and fossil groups, and see if these systems share a common 
origin. 

In this p aper, we use the Millennium Simulation 



Springel et al.ll 2005) together with a semi-analytical model 



Lucia fe Blaizot 20071 ) of galaxy formation within dark 



matter haloes to identify isolated field elliptical galaxies, to 
study their properties and formation history, and to com- 
pare them with observational data. This paper is organised 
as follows. In Section [21 we discuss the Millennium Simula- 
tion, the semi-analytical galaxy catalogue used, and sample 
selection. In Section|3j we compare the properties of the Mil- 
lennium IfEs with observations, while in Section [3] we show 
that the IfEs form a population that is different from the reg- 
ular ellipticals. In Section [5] we concentrate on evolution of 
IfEs and possible formation mechanisms, while in Section [5] 
we discuss our results, the formation of field elliptical galax- 
ies and their possible connection to fossil groups. Finally, 
in Section we summarise our results and draw the con- 
clusions. Throughout this paper, we adopt a parametrized 
Hubble constant: H = WOh km s" 1 Mpc -1 . 



2 SAMPLE SELECTION 

2.1 The Millennium Simulation 

We use a simulation that covers a big enough spatial vol- 
ume and has a sufficient mass resolu tion - the dark mat - 
ter only Millennium Simulation (MS: ISpringel et al]|2005l ). 
a 2160 3 -particle model of a co-moving cube of size 500/i _1 
Mpc, on top of which a publicly av ailable semi-analytical 
galaxy model (|Lucia fc B laizot 20(53) has been constructed. 
The cosmological parameters used in the MS were: Q m — 
ttdm + n b = 0.25, Q b = 0.045, h = 0.73, Q A = 0.75, n = 1, 
and us = 0.9 where the Hubble constant is characterized as 
100ft km s -1 Mpc" 1 (for detailed description of the MS, see 
ISpringel et alj |2005j). These values were inferred from the 
first-year W MAP (Wilkinson Mic rowave Anisotropy Probe) 
observations l|Spergel et al]|2003l ). 

The galaxy and dark matter halo formation modeling 
of the MS data is based on merger trees, built from 64 
individual snapshots. From these time-slices, merger trees 
are built by combining the information of all dark matter 
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haloes found at any given output time. This enables us to 
trace the formation history and growth of haloes and sub- 
haloes through time. Properties of galaxies in the MS data 
are obtained by using semi-analytic galaxy formation mod- 
els (SAMs), where star formation and its regulation by feed- 
back processes is parametrized in terms of analytical phys- 
ical models. The semi-analytical galaxy catalogue we use 
contains about nine million galaxies at z = down to a 
limiting absolute magnitude of Mr — 51og/i = —16.6. A 
detailed description of the creation of the MS Galaxy Cat- 
alogue , used in this study, can be found in lLucia fc Blaizot] 
l|2007l ). 

Semi-analytical galaxy formation models are known to 
be less than perfect. There are several free parameters in 
each SAM that all have a different impact on the properties 
of galaxies. The free parameters usually control the feedback 
effects and cooling of hot gas. The SAM also controls how 
gas is stripped from a galaxy when it approaches another 
galaxy, and possible starburst events. The physics of these 
processes described is not well known. However, in general, 
SAMs can reproduce the observed galaxy luminosity func- 
tion and other statistics well. As our study is statistical by 
nature, it is likely that the SAM used for the Millennium 
Simulation Galaxy Catalogue is accurate enough to predict 
the properties of IfEs and their evolution and formation his- 
tory. 

The MS Galaxy Catalogue does not predict galaxy mor- 
phologies. To assign a mo rphology for every galaxy, we use 
their bulge-to-disk ratios. ISimien fc de Vaucouleursl (|l986h 
found a correlation between the B— band bulge-to-disc ratio 
and the Hubble type T of galaxies from observations. The 
mean relation may be written: 

< Am(T) >= 0.324x(T) - 0.054x(T) 2 + 0.0047x(T) 3 , (1) 

where Am(T) is the difference between the bulge magnitude 
and the total magnitude and x(T) = T + 5. We classify 
galaxies with T < —2.5 as ellipticals, those with —2.5 < 
T < 0.92 as SOs, and those with T > 0.92 as spirals and 
irregulars. Galaxies without any bulge are classified as type 
T = 9. 



2.2 Selection of IfEs 

There are a few different selection criteria that are being 
used to define isolated field elliptical galaxies in observations 
(seee g.lColbert et al.l200ll ; lMarcum et alj20oH ; lReda etafl 
12004 ICollobert et alj|2006h . Despite the differences in selec- 
tion criteria, all studies identifying IfEs share a common 
ideology. In general, IfEs are defined as elliptical galaxies 
that do not have nearby optically bright companion galax- 
ies, usually in the B— band. This is tested by using an iso- 
lation sphere or a cone and choosing a minimum magnitude 
difference between the brightest and the second brightest 
galaxy Am^ inside the sphere or the cone. In observational 
studies, the isolation cone is often taken as a circle in the 
sky that expands in redshift space. Due to peculiar veloci- 
ties, accurate line-of-sight distances of galaxies are unknown 
complicating the identification. Thus, observational studies 
differ in the numerical values of the isolation criteria. De- 
spite these differences, discussed next, we use all possible 
available observational data for comparison as the number 
of observed targets is small. For completeness, we show in 



Section T3. II how the number of IfEs depends on the selection 

crit eria. 

IColbert et al.l (|200lf ) adopted rather strict rules for iso- 
lation, as they required that IfEs cannot have catalogued 
galaxies with known redshifts within a projected radius of 
1.0/i _1 Mpc and a velocity of ±1000 km s _1 , while for the 
morphology of the IfEs they required Hubble types T ^ —3. 
However, as their source of galaxies w as the Third Reference 
Cata logue of Bright Galaxies (RC3, Ide Vaucouleurs et al.l 
1991), their sample IfEs might have faint companion galaxies 
due to the incompleteness of the source catalogue. Missing 
companions do not affect the study of bright isolated galax- 
ies per se, but their environment is affected by the incom- 
pleteness and missing companions. iMarcum et alj (|2004l ) 
based their initial sample of IfEs on the same catalogue. 
They adopted an even stricter criterion in projected radius, 
a minimum projected physical distance of 2.5 Mpc to any 
neares t neighbour brighter than My = —16.5. iReda et al.l 
(2004) used less strict criteria in their study and required 
only that an IfE candidate has no neighbours with bright- 
ness difference Amu ^ 2.0 mag, in the B— band, within 
0.67 Mpc in the plane of sky and 700 km s _1 in recession 
velocity. Further, their elliptical galaxies are also of Hubble 
type T < -3. 

For the selection of IfEs fro m the Millenn i um S imulation 
we adopt the criteri a used by Smith et al.1 |2004l 'l. similar 
to those adopted by IReda et al. I (|2004h .We concentrate on 
galaxies that are incontrovertibly ellipticals and thus adopt 
a strict morphology limit T ^ —4. We further limit the 
brightness of the IfE candidates with the magnitude cut- 
off of Mb ^ —19. For the isolation criterion we adopt two 
isolation spheres with the radii of 0.5/i -1 and 1.0ft _1 Mpc. 
Within these isolation spheres, we require that the B— band 
magnitude difference of the first- and second-ranked galax- 
ies Amu must be 2.2 and ^ 0.7 ma g, respectively. The 
application of the criteria is illustrated in lSmith et al.l (|2004l . 
Fig.l). 

Our adopted values of magnitude differences correspond 
to factors of about eight and two in luminosity for the small 
and large sphere, respectively. This choice ensures that pos- 
sible companions are small and light enough not to produce 
any major perturbations in the gravitational potential of 
the system. Note, that in case of simulations our isolation 
criteria operate in real space rather than in redshift space. 
Therefore, we can use isolation spheres rather than cones 
and we are not plagued by interlopers due to inaccuracies in 
distance measurements. 



2.3 Simulated IfEs and Control Samples 

To identify isolated field elliptical candidates in the MS, we 
chose five independent cubic volumes inside the simulation 
box. Initially, we chose each volume to have a side length of 
~ 200/i.- 1 Mpc, while none of the cubes overlap. However, 
we later limited the volume from where the possible IfE 
candidates can be identified to have a side length of 195/i _1 
Mpc. This was done to enable the study of the surroundings 
of each IfE candidate in the same fashion. Moreover, without 
the limitation, we could have accidentally identified a field 
elliptical candidate ^ 1. Oft -1 Mpc from an edge of the box, 
leading to a possible false identification. The five volumes 
were chosen to overcome computational issues and for easier 
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study of IfE environment. With five independent volumes it 
is also possible to make comparisons between the IfEs of 
each volume. 

At first all galaxies inside each volume are treated as 
possible IfE candidates. After identifying all candidates with 
T ^ —4 and Mb ^ —19 we apply the isolation criteria 
discussed in the previous Section. This produces an initial 
list of 302 galaxies that fulfill the criteria for isolated field 
elliptical galaxies. After compiling the initial list of IfEs, 
we inspected every galaxy individually, and found that nine 
candidates were not the main galaxy of the dark matter 
halo they reside in. Therefore, all galaxies that fulfill the 
observationally motivated IfE criteria are not the central 
galaxies of their Friends-of-Friends dark matter groups. In 
a strict sense these are not isolated galaxies, as they belong 
to a dark matter halo containing galaxies more massive and 
more luminous than the IfE candidate. Therefore, we omit 
them from the final list of IfEs, which contains 293 galaxies. 

We also compile two control samples of elliptical galax- 
ies for comparisons. The first control sample contains in total 
4563 elliptical (T sC —4) galaxies brighter than —19 in the 
B— band at the redshift z — and it is named as Ellipticals 
(abridged as Es). For the second control sample, called Main 
Ellipticals (abridged as MEs) , we only select galaxies of the 
first control sample that are the main galaxies of their dark 
matter haloes. This requirement further limits the Ellipti- 
cals sample, and we are left with 1209 galaxies in the Main 
Ellipticals sample. 

2.3.1 Non-standard IfEs 

The removed galaxies are interesting from another point 
of view than regular IfEs, as these nine non-standard IfEs 
are 'subhalo' galaxies (i.e. satellites of a larger galaxy) and 
have multiple nearby companions, from 30 to 80 inside the 
large isolation sphere. Thus, these galaxies are not isolated 
in galaxy number density, but reside in a cluster rather than 
in the field. Physically it is obvious that a satellite galaxy is 
not the dominant galaxy in its dark matter halo. Since we 
find these non-standard IfEs, our observationally motivated 
isolation criteria are not strict enough and the radii of the 
isolation spheres should be slightly larger to avoid any false 
identifications. But as the number of non-standard IfEs is 
small, we can keep the criteria and remove the non-standard 
galaxies from our final sample, as described above. 

The satellite non-standard "isolated" ellip tical galaxies 
can s hed some light on the observational result (|Smith et al.l 
|200ST ) where an IfE galaxy (NGC 1600) was also found to 
be located of f from the dynamical centre of the system. 
ISivakoff et ail (|2004h found that the X -ray emission is cen- 
tered slightly to the northeast of NGC 1600, suggesting that 
the galaxy is not at the centre of the gravitational poten- 
tial. The dynamical study and the X-ray observations to- 
gether suggest that NGC 1600 is surrounded by a massive 
halo extending out to several hundred kiloparsecs. A sub- 
halo galaxy, fulfilling the optical IfE criteria, could explain 
the observations of NGC 1600. All nine galaxies in our sim- 
ulations reside inside a large dark matter halo, and in ob- 
servations would be found to be off from the centre of the 
potential well, as in the case of NGC 1600. However, a mas- 
sive and more luminous galaxy than NGC 1600 has not been 
detected at the centre of the potential well in observations, 



while this is often the case in simulations. Thus, a more 
straightforward explanation, where the IfE orbits around 
the centre of the potential well and therefore seems to be 
shifted from the centre, seen in X-ray observations, is more 
plausible. This is in agreement with an X-ray study of an- 
other isolated elliptical NGC 4555, which is assumed not to 
lie in the centre of a massiv e group-scale dark matter halo 
l|Q'Sullivan fc Ponmanlliobl ). 



3 PROPERTIES OF OBSERVED AND 
SIMULATED FIELD ELLIPTICAL 
GALAXIES 

Below we compare the properties of simulated IfEs with 
those of the observed IfEs. If these are close enough, it will 
justify the theoretical study of the formation, evolution and 
merger histories of IfEs, using the simulated IfE population. 
Due to the large number of simulated IfEs we can actually 
make theoretical predictions for properties, formation mech- 
anisms and evolution of IfEs that can be tested against ob- 
servational data when the sample of observations is large 
enough. 

In following sections we use the Kolmogorov— Smirnov 
(KS) two-sample test to study the maximum deviation be- 
tween the cumulative distributions of two samples. The null 
hypothesis of the KS test is that the two samples are from 
the same population. We give our results as probabilities (p- 
values) that the difference between the two samples could 
have arisen by change if they are drawn from the same par- 
ent population. 

3.1 Number of IfEs 

The fraction of isolated field elliptical (T ^ —4 and Mb ^ 
— 19 mag) galaxies among all galaxies of any brightness, mor- 
phology or dark matter halo status in the MS is merely 
~ 3.5 x 10 -3 per cent, corresponding to a number density of 
~ 8.0 x 10 _6 h 3 Mpc -3 . This is a very low number density, 
however, such a straightforward calculation is not totally 
justified. More meaningful number to conside is the fraction 
of IfEs among elliptical galaxies. Moreover, our morphology 
limit (T ^ —4) for a simulated IfE is very strict, biasing the 
fraction of elliptical galaxies to significantly lower fractions 
than observed in the real Universe. Thus, we also identi- 
fied isolated field elliptical galaxies with relaxed morphology 
limit T < —2.5 and quote the fractions below. 

About 0.19 per cent of all simulated elliptical (now T < 
—2.5) galaxies of any brightness, can be identified as IfEs, 
when the criteria of Section 12.21 are adopted for IfEs. If we 
relax the strict morphology limit of the IfEs and require that 
IfEs also have T < —2.5 the fraction rises to ~ 0.48 per cent 
as about 2.5 times more IfEs can be identified. Relaxing the 
morphology limit will therefore change the number density 
of IfEs to ~ 1.9 x 10" 5 h 3 Mpc" 3 . When the brightness of 
the simulated elliptical galaxies and their morphologies are 
limited to Mb ^ —19 at z = and T ^ —4, respectively, 
about 6.4 per cent of these galaxies are identified as IfEs 
at redshift zero. If we further limit our simulated elliptical 
galaxies to galaxies that are the main galaxies of their dark 
matter haloes (as in case of IfEs), ~ 32 per cent of ellipticals 
are now identified as IfEs. Thus, IfEs are very rare objects 
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among all galaxies, but at the same time, over 30 per cent of 
main elliptical galaxies brighter than —19 mag in B— band 
can be classified as IfEs. 

IStocke et all |2004 ) found from magnitude limited 
(ms ^ 15.7) observations that less than three per cent of all 
galaxies can be classified as isolated. This value should be 
interpreted as an upper limit for t he number o f IfEs, as it 
was based on all types of galaxies. IStocke et al. | l|2004) fur- 
ther found that early-type galaxies outnumber SO galaxies 
2 : 1 in very isolated areas. However, as spiral gal axies are 
the m ost dominant galaxy type in low densities fe.g. lDressler] 
1980), the observable estimate of t he number of IfEs is sig- 
nifican tly lower than one per cent. iHernandez- Toledo et al.l 
(2008) concluded that early-type galaxies amount to only 3.5 
per cent among all isolated galaxies, lowering the observa- 
tional estimate of IfEs to ~ 1 x 10~ 2 per cent. This fraction 
is in good agreement with our findings as about ~ 9 x 10~ 3 
per cent of all simulated galaxies can be identified as IfEs 
when T < —2.5 morphology limit has been adopted. 

As observational studies adopt different criteria for iso- 
lation, we illustrate their influence on the number of IfEs, 
with a stricter set of parameter values. For simplicity, we 
adopt only one isolation sphere with a radius of 2.5h~ 
Mpc, and require that the magnitude difference between the 
brightest and the second brightest galaxy Ami2 is ^ 2.0 mag 
inside the isolation sphere. For the Hubble type, we adopt 
the same requirement as previ ously for our I f Es: T ^ —4. 
Such values were used, e.g., bv lMarcum et all (|2004T ). With 
these values we find that only ~ 3 x 10 -4 per cent of all MS 
galaxies can be classified as IfEs, corresponding to an ex- 
tremely low number density of ~ 6 x 10~ 7 h 3 Mpc -3 . This is 
approximately ten times less than before. It is obvious that 
changing the isolation criteria has a major impact on the 
number of IfEs and therefore their statistical properties, as 
stricter parameter values require a galaxy that is located in 
a true void. This result also shows that galaxies with compa- 
rable luminosities tend to reside in relatively close proximity 
rather than being truly isolated from any other reasonable 
sized companion. 

The results of this section show that the fraction of iso- 
lated field elliptical galaxies depends highly on the isolation 
criteria adopted, but also on the morphological type that is 
required. As the morphological type of simulated galaxies is 
less than precise our number density for IfEs should not be 
taken as a strict limit, but as a guideline when observations 
are being planned. Thus, higher fractions of IfEs can be ex- 
pected if isolation criteria or morphology limits are being 
relaxed. 

3.2 Colour-magnitude Diagrams 

The easiest property of IfEs to observe must be the lumi- 
nosity, and thus the colour-magnitude relation. This simple 
yet powerful relation is studied in Figs. [T] and [2] which show 
the colour-magnitude diagrams (CMDs) of simulated and 
observed IfEs. Fig. [1] shows a comparison between the sim- 
ul ated IfEs and obse rved field elliptical galaxies presented 
in iReda et al.l (|2004T ). while Fig . [2] shows a comparis on to 
the sample of IfEs presented in iMarcum et al.l (|2004l 'l . We 
also plot the Ellipticals (our control sample) for complete- 
ness. Note, however, that in this section we limit the discus- 
sion of the colour-magnitude diagrams to the comparison of 



Table 1. Colours of simulated and observed IfEs and simulated 
non-isolated elliptical galaxies (Es). 



Sample 


Amag 


Mean Colour 


Error 


17 


IfEs 


B - 


R 


1.47 


0.01 


0.23 


IfEs 1 


B - 


R 


1.57 


0.01 


0.17 


IfEs 


B - 


V 


0.79 


0.01 


0.15 


IfEs 2 


B - 


V 


0.82 


0.01 


0.14 


Es 


B - 


R 


1.58 


0.002 


0.10 


Es 


B - 


V 


0.86 


0.01 


0.07 


Reda et al. (2004) 


B - 


R 


1.49 


0.06 


0.16 


Marcum et al. ("2004) 


B - 


V 


0.76 


0.06 


0.18 



Note: The error refers to the standard error of the mean while 



a is the standard deviation. IfEs 1 and IfEs 2 refers to samples 
where faint simulated isolated field elliptical galaxies have 
been removed and the IfE samples contain only galaxies with 
Mr < -21.5 and M B < -19.5, respectively. 



simulated and observed IfEs, while the differences between 
simulated IfEs and Es are discussed in Section \4. 21 Colours 
and associated errors and standard deviations of observed 
and simulated IfEs and simulated elliptical galaxies (Es) are 
listed in Table □ 

Fig. [T] shows that the observed IfEs lie in slightly bluer 
region in the CMD than simulated IfEs of the same bright- 
ness and that they seem to follow a tighter correlation 
than simulated IfEs. The B — R colour scatter of simu- 
lated IfEs is slightly larger compared to observed IfEs: all 
simulated IfEs sho w a la dev i ation of 0.23, however, as 
the observations of IReda et all (|2004l ) are limited to rela- 
tively bright (Mr ^ —21.5) galaxies, we calculate the colour 
scatter of bright simulated IfEs using magnitude cut-off of 
Mr ^ —21.5 and find it to be 0.17. The colour scatter of the 
bright simulated IfEs is close to the observed B — R colour 
scatter (0.16), yet slightly higher. However, one should bear 
in mind that the num ber of observed IfEs is still small. 

IReda et~afl (|2004r ) found a mean effective colour of (B — 
R) e — 1.49 ±0.06 for their isolated elliptical galaxies, which 
is in good agreement with our mean value of 1.47±0.01 (the 
errors are the standard error on the mean, see Table [1] for 
details). If we limit our sample of simulated IfEs to galaxies 
bright er than Mr ^ —21.5, as in the sample of Red a et al.l 
(2004), the mean B - R colour be comes 1.57 ± 0.01. This 
is redder than the mean colour in IReda et al l l|2004l ) and 
shows that our mean colour is significantly affected by faint 
IfEs of th e blue cloud. Whe n the KS test is applied to the 
observed (iReda et al.l I2004T ) and all simulated IfE colours 
the p-value (~ 0.33) shows that we cannot reject the null 
hypothesis at 30 per cent level. However, if we consider only 
bright (Mr ^ —21.5) simulated IfEs, the p-value is only 
~ 0.01, implying that a difference between the observed and 
simulated IfE B — R colours exists when faint simulated IfEs 
are excluded from the comparison. 

A straigth line fit to both observed and simulated IfEs 
differs in the slopes and intercepts (see Fig. [TJ . The linear fit 
of observed IfEs shows a steeper slope tha n the fit of the sim - 
ulated IfEs. However, as the IfE sample of I Reda et all (|2004f i 
miss galaxies fainter than Mr > —21.5 it is not entirely clear 
what the order of magnitude of this difference might be, and 
due to the small number of observed IfEs in the sample of 



6 Sami-Matias Niemi et al. 



iReda et al.l l|2004h only a single new data point at the faint 
end of the CMD could change the fit significantly. 

The V-band (Fig. CMD differs from the 7?-band 
diagram as correlations seem significantly looser. Here the 
colour scatter of the observed IfEs is larger than for the sim- 
ulated IfEs (la de viations are 0.18 and 0.15, respectively). 
The IfE sample of iMarcum et alJ (|2004l ) extends to slightl y 
fainter magnitudes than the sample of IReda et al.l {2004), 
but their fainter IfEs are surprisingly red (B — V ~ 0.9). The 
simulated IfEs show slightly redder colours than observed 
IfEs; the mean B — V colours are 0.79 ±0.01 and 0.76 ±0.06, 
for simulated and observed IfEs, respectively. The colours 
of simulated IfEs agree with observed colours within their 
standard errors of the mean, and the KS test approves the 
null hypothesis with high probability (p-value ~ 0.63) when 
colours of the simul ated IfEs are comp are to the colours of 
the IfE sample of IMarcum et alJ | |2004). The p-value would 
be only 0.04 if the observed IfEs of IMarcum et al.l (|2004l ) 
were compared against the simulated Es. T hus, the simu- 
lated IfEs agree well with the observed IfEs of lMarcum et al.l 
(|2004h when B — V colour samples are compared. 

Figs. [T] and [2] show a separate population of simulated 
IfEs that belong to the blue cloud and populate the faint end 
of the CMD. Quantitatively the colour and the brightness of 
this population is following: B — R ^ 1A or B — V ^ 0.7 and 
Mr > -21.5 or Mb > -20.0. The population of faint and 
blue IfEs comprise ~ 26 per cent (76) of all simulated IfEs, 
thus, every fourth IfE belong to this po pul ation. We note 
that n one of the IfEs of IReda et~ail (|2004l ) or IMarcum et all 
(2004) po pulates this fain t and blue part of the CMDs. The 
sample of Mar cum et al.l (|2004l ) shows two very blue IfEs 
that are part of the global blue cloud, but are not part of 
the population of faint and blue IfEs that is found from the 
simulati on. However, this mi ght be due to the magnitude 
limit of lMarcum et all l|2004h sample, as they do not have 
IfEs fainter than M B > - 19.5. 

Im 

arcum et a l. (2004) found preliminary evidence that 
50 per cent of their sample of isolated early-type galaxies 
show blue global colours. However, because of the small sam- 
ple size they could not conclude that the higher occurrence 
of blue E-type galaxies is related to the extremely low den- 
sities of the associated environments. As we are not limited 
by a small sample, we can confirm whether IfEs show bluer 
global colours than non-isolated early-type galaxies or not. 
The mean B — V colour of sim ulated IfEs is ~ 0.79 ± 0.01 in 
agreement with the findings of lMarcum et al' (2004). How- 
ever, this value is significantly affected by the population of 
faint IfEs; if we remove the faint IfEs and consider only IfEs 
with Mb ^ —19.5 the mean colour changes to ~ 0.82 ±0.01 
(see Table [T] and also Fig. [6] of the next Section) . This is 
close, yet slightly bluer, than the mean B — V colour of sim- 
ulated ellipticals (0.86 ± 0.01). A similar result is seen when 
the B — R colours are studied; bright (Mr ^ —21.5) sim- 
ulated IfEs are almost as red as all simulated non-isolated 
ellipticals (Es) B - R ~ 1.57 ± 0.01 and 1.58 ± 0.002, re- 
spectively. Thus, our results show that simulated IfEs do 
not show global blue colours if faint IfEs are removed, but 
the blue colours are due to the separate, faint and blue pop- 
ulation of IfEs. Our theoretical findings predict that ~ 26 
per cent of IfEs show global blue colours and that these IfEs 
belong to the separate population of faint and blue IfEs. 
Different trends in the CMDs, visible in Figs. [1] and [2l 
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Figure 1. Colour (B — R) versus absolute R— band magnitude 
for the observed and simulated IfEs. The data for the Es control 
sample have been plotted for comparison. For IfEs and Es a linear 
fit is shown. For clarity, only every 10th Es has been plotted. Note 
the separate population of faint and blue IfEs. 
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Figure 2. Colour (B — V) versus absolute B— band magnitude 
for the observed and simulated IfEs. The data for the Es control 
sample have been plotted for comparison. For Es a linear fit is 
shown. For clarity, only every 10th Es has been plotted. Note the 
separate population of faint and blue IfEs. 



may be interpreted as environment effects of galaxy forma- 
tion. However, it is also possible that the bluest IfEs have 
had different formation mechanism and evolutionary path 
than redder IfEs. Different evolutionary paths and merging 
histories could explain different properties observed at red- 
shift z — 0, as well as the influence of environment. Both 
possibilities will be explored in detail later. 



3.3 Dark Matter Halo Masses 

We continue comparisons of simulated and observed IfEs by 
studying their dark matter haloes. Even it is far from simple 
task to derive reliable dark matter halo masses from observa- 
tions the comparison is highly interesting, as the dark matter 
halo properties and galaxy properties are tightly linked. 
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Isolated field elliptical galaxies in the MS are mainly 
found residing inside dark matter haloes that are lighter 
than 7 x 10 12 h~ 1 MQ. The dark matter haloes of simu- 
lated IfEs are surprisingly light; the median mass is only 
~ 1.2 x 1O 12 /i -1 M0. Even the most massive dark mat- 
ter halo hosting an IfE is lighter than 2.2 x lO 13 fr _1 M , 
which is comp a rable to a dark matter halo of a small group. 
iMemola et at! (|2009l ) calculated the total masses of two of 
their isolated ellipticals NGC 7052 and NGC 7785 from 
X-ray observations and quote values ~ 5 x lO^M© and 
~ 1.9 x 1O 12 M0, respectively. These mass va lues agree well 
with o ur findings of dark matter halo mass. iNorberg et al.l 
(2008) did not find any isolated systems residing in haloes 
outside the range ~5x 10 11 to lO 13 /i _1 M0 in their simula- 
tions, therefore, we can conclude that IfEs reside in lighter 
than ~ 2 x 10 13 /i _1 M Q dark matter haloes. 

The dark matter halo mass also sets constrains to the 
possible formation mechanisms; IfEs residing in light haloes 
cannot be merger remnants of groups or clusters. Only those 
IfEs that have a large and massive dark matter halo could 
have formed via merger of a group or multiple larger galax- 
ies, but even then the group should have been poor with 
only few member galaxies. Evidently, this is true for fos- 
sil groups which are in many ways similar objects to IfEs. 
The dynamical mass es of fossil groups range fro m ~ 10 13 to 
lO 14 h -1 M (see e.g. iDi'az-Gimenez et ai]|200Sf ). 

3.4 Ages 

The last property of IfEs we compare is their age. The me- 
dian mass weighted age of our model IfEs is 8.84 Gyr (see 
Table [3]). We also note that the number density of young 
(mass weighted age less than 5 Gyr) IfEs is extremely low 
(Fig. 1101) . g iving a low er limit for the age of an IfE. 

iReda etal] (|2005l ) found a mean age o f their IfEs to 
be 4.6 ± 1.4 Gyr, while iProctor et all (2005) quote an age 
estimate of ~ 4 Gyr. These age estimates are around the 
young end of our values. However, these estimates should 
not be compared directly to our values as the definitions 
are different: mass- vs . luminosity-weighted age. Moreover, 
ICollobert et"afl l|2006l ) found a broad range of stellar ages 
for their IfEs; ranging from ~ 2 to 15 Gyr, in good agree- 
ment with our estimates, as the mass weighted age of model 
IfEs ranges from 5 to 12 Gyrs. The big scatter suggests that 
the formation and evolution of IfEs is not concentrated at a 
fixed epoch. This can also explain the large scatter we see 
in some properties of IfEs, as different formation times and 
evolutionary paths can lead to significantly different prop- 
erties at z — 0. 



3.5 Environments 

To get a complete picture of the properties of isolated field 
elliptical galaxies, we have to look at their environment. We 
define a companion as a galaxy that resides inside either of 
the two isolation spheres; the small 0.5/i -1 Mpc or the large 
sphere 1.0ft -1 Mpc. If not specifically mentioned, the num- 
ber of companions refers to the total number of galaxies, 
Ncomp, inside the large isolation sphere. Our definition of 
a companion does not guarantee that the galaxy belongs to 
the same dark matter halo as the IfE. Thus, it is possible 



that some of the companion galaxies are not gravitationally 
bound to the IfE. Truly isolated galaxies, without any com- 
panions in simulations, should be treated with caution. It 
may be that the mass resolution is not sufficient to form 
small subhaloes that could harbour a dwarf galaxy. 

IfEs can have from to ~ 30 companion galaxies. We 
do not find a single IfE with more than 30 dwarf compan- 
ions, indicating that our IfEs are in relatively low density 
environments. The majority of IfEs have to 20 compan- 
ions, while the mean value of the number of companions, 
Ncomp, is 10.7 in the MS. For IfEs with the total number 
of companion galaxies less than 17, most of the companions 
are found to reside inside the smaller (0.5ft _1 Mpc) isola- 
tion sphere. As the magnitude difference limit inside the 
small sphere is 2.2 we can be sure that our IfEs (excluding 
the nine subhalo galaxies discussed earlier) are well isolated 
from bright nearby galaxies. So, IfEs do not reside inside 
cluster-sized dark matter haloes, with virial radii ~ 1.5/i -1 
Mpc, as most their companion galaxies are within 0.5ft~ 
Mpc distance from the halo's main galaxy. If IfEs resided 
in cluster haloes, we should find companion galaxies also in 
the larger isolation sphere. We would also expect to find a 
larger number of dwarf companions. 

IReda et al. found that only the very faint dwarf 

galaxies (Mr ^ —15.5) appear to be associated with iso- 
lated ellipticals. On the contrary, we find companion galax- 
ies with a broad range of magnitudes from the MS. The 
mean R— band magnitude for these companions is —17 mag. 
The quartile values for companion R— band magnitudes are 
— 18.2 and —15.8 mag showing that model IfEs can have rel- 
atively bright companion galaxies. Thus, the magnitude gap 
between companion galaxies and an IfE seems to be larger 
in observations than in simulations. This may be due to the 
way how galaxy luminosities are treated in the SAM. Thus, 
this result is probably not without a bias due to the limiting 
mass resolution of the MS. 

In Fig. [3] we show the distribution of companion galaxy 
distances from IfEs. We have divided the IfEs sample into 
two, to separate the population of blue, light and faint IfEs; 
the IfE sample is divided by dark matter halo mass. Fig. 
[3] clearly shows that more massive IfEs have close compan- 
ions more often than light IfEs. The mean virial radius of 
the dark matter haloes of heavy IfEs (Mdm > 1O 12 M0) is 
225ft" 1 kpc, while it is 115ft" 1 kpc for light IfEs (Mdm < 
1O 12 M0). A significant number of companions of heavy IfEs 
are inside the mean virial radius and belong to the same 
dark matter halo as the IfE. However, for blue, light and 
faint IfEs the trend is opposite; we find most of the com- 
panions more than six times the mean virial distance away 
from the IfE. The median companion distances are 0.33 and 
0.73ft _ Mpc for heavy and light IfEs, respectively. 

Similar results are also found if the IfEs are divided 
by their colour. The blue, B — R 1.4, IfEs have most of 
their companion galaxies more than five times away than 
the mean virial radius, while the red IfEs have the majority 
of their companions within the mean virial radius. These 
results show that the most isolated field elliptical galaxies 
are blue and relatively faint. 
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Figure 3. Total number of companion galaxies, Nqomp, inside 
the large isolation sphere for light and massive IfEs. 



4 COMPARISON OF REGULAR AND FIELD 
ELLIPTICALS 

Here we compare the properties of the simulated isolated 
field elliptical galaxies to the control sample of all elliptical 
galaxies in the simulation, to find out the differences between 
the two populations. For many properties the median val- 
ues are relatively close to each other, but clear differences 
exist in the shapes of the distributions. Thus, throughout 
the following sections we use the Kolmogorov— Smirnov two- 
sample test to assess whether two distributions are drawn 
from the same parent population. Results of the KS tests 
are given as probabilities (p-values) and presented in Table 
[2] For completeness, we also plot the histograms of studied 
properties (see Figs. [5] — 110[) and present the quartile values 
in Table El 



4.1 Morphology 

The distributions of morphologies (T values) of isolated field 
elliptical galaxies and Es are very similar, with median val- 
ues: —5.62 and —5.80, respectively. The KS test approves 
the null hypothesis with high probability (p- value ~ 0.39), 
thus, it is likely that the differences in morphology distri- 
butions have arisen by change. This is no surprise as our 
definition of IfEs and Es requires that the morphology value 
is smaller than four. 



4.2 Galaxy Colour-magnitude Diagrams 

The galaxy colour-magnitude diagrams show three main fea- 
tures: the red sequence, the green valley, and the blue cloud. 
In general, elliptical galaxies populate the area known as the 
red sequence. This is certainly true for most elliptical galax- 
ies of the MS, however, is this true for IfEs as well? 

Figs. Q] and [2] show that the red sequence of IfEs starts 
at Mr ~ —21 (or Mb ~ —19.5), while it extends to fainter 
magnitudes for Es. A trend of redder colours with increasing 
luminosity is noted for both simulated IfEs and Es. However, 
the slope of the trend is steeper and the scatter is higher for 
IfEs. Simulated IfEs show a broader distribution in colours 



(Figs. [T] and [2]) than simulated Es. The la scatter of B — R 
colours is ~ 0.23 and ~ 0.10 for IfEs and Es, respectively. 
For B — V colours the la values are ~ 0.15 and ~ 0.07 for 
IfEs and Es, respectively. Spearman rank order correlation 
test shows a very strong (correlation coefficient cc ~ —0.65) 
correlation for IfEs (in Fig. [TJ and significant (cc ~ —0.33) 
correlation for Es. The trends in Fig. [2] are not as strong 
according to the Spearman test: cc ~ —0.60 and ~ —0.27 for 
IfEs and Es, respectively. However, all correlations are highly 
significant, as the probability to have as large correlation 
coefficients for uncorrelated data is < 10~ 20 in all cases. 

The colour correlations in Figs. Q] and [5] suggest that 
IfEs follow a different colour-magnitude trend than Es. The 
KS tests show large differences (D-values) when the colours 
of IfEs are compared to Es. The KS test rejects the null 
hypothesis in case of IfEs and Es B — R and B — V colours 
with high probability (p-values < 10~ 5 in both cases). If we 
exclude the separate population of faint and blue IfEs, the 
previous statement does not change, only the probabilities 
(p-value now < 10~ B ). Thus, the CMDs of IfEs and Es differ 
significantly even when the separate population of faint and 
blue IfEs is removed. 



4.3 Colour-mass Diagrams 

Fig. U shows the B — R colour of a galaxy as a function of 
the mass of the underlying dark-matter halo. We note from 
Fig. [4] that the Es show a rather constant relation with a 
few outliers; in general, Es have the B — R colour ~ 1.6 ± 
0.15 independent of the mass of the dark matter halo they 
reside in. However, a completely different trend is seen for 
IfEs, as they tend to get redder when the dark matter halo 
mass grows, while light (Mum < 10 12 ft _1 M©) dark matter 
haloes host blue IfEs. The B — R colours of IfEs grow steeply 
as a function of dark matter halo mass below 1O 12 /i _1 M 
indicating the influence of dark matter halo properties to 
the IfE galaxy they host. 

All IfEs are the main galaxies of their dark matter 
haloes. If instead of the Es sample we used the MEs for 
comparison, Fig. [4] would change. There would still be a 
significant number of red, B — R ~ 1.6, non-isolated el- 
liptical galaxies residing in dark matter haloes lighter than 
Mdm < 10 12 /i _1 Mq, but their total number would be sig- 
nificantly smaller. MEs that reside in light {Mdm < lO 12 /!" 1 
Mq) dark matter haloes show a big scatter in colours, in- 
dicating that not only the dark matter halo, but also the 
environment and the formation history of the halo effects 
the galaxy it hosts. 

Fig. [4] further shows that the population of faint and 
blue IfEs noted in the colour-magnitude diagrams resides in 
the lightest dark matter haloes with Mdm < 10 12 /i _1 Mq. 
Below this mass scale we see a strong correlation between 
the dark matter halo mass and the evolution of the galaxy. 
Thus, simulations predict an unobserved population of IfEs 
that are blue and faint and that reside in light dark matter 
haloes. 



4.4 Colour and Luminosity Distributions 

Fig. [5] shows the number density of absolute B— band rest 
frame magnitudes for IfEs and Es. The distributions look 
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Figure 4. Dark matter halo mass versus colour (B — R) of the 
galaxy. For clarity, only every 5th Es has been plotted. 



Figure 6. Number density of colour distribution (B — R magni- 
tudes) for isolated field ellipticals (IfEs) and Es. 
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Figure 5. Number density of B— band rest frame magnitude for 
isolated field ellipticals (IfEs) and Es. 

very different and the KS test (Table [2} rejects the null 
hypothesis with high probability. The differences in the lu- 
minosity functions are likely due to the isolation criteria 
of IfEs. The most interesting result, seen in the Figure, is 
that isolated field elliptical galaxies have an almost constant, 
~ 8 x 10~ 6 , number density throughout their B— band mag- 
nitudes. Surprisingly, when only brighter (—21.7 < Ms < 
— 21.0) elliptical galaxies are considered, it is almost as prob- 
able to find an isolated field elliptical as it is to find a non- 
isolated elliptical. It is also noteworthy that we have not 
identified a single IfE galaxy brighter than —21.7 mag (in 
B— band), while the brightest Es are almost one magnitude 
brighter (~ —22.5). This result is natural, since the bright- 
est elliptical galaxies are normally found to reside in centres 
of large clusters. The brightest cluster galaxies (BGCs) have 
My ~ —23.5 mag and they have usually e xperienced many 
merg ing events during their evolution (e.g. iLucia &: Blaizotj 
120071 ). Multiple merging events have brought more mass and 
gas to the central galaxy and have caused massive star for- 
mation and greater luminosity. It is unlikely that IfEs have 
followed the same evolutionary path (see Fig. [7| and [8]) . 



Since the B— band luminosity function of IfEs is almost 
constant, it leads to another question, namely if IfEs have a 
constant number density in colour as well. Fig. [6] replies to 
this question and presents the number density for the B— R 
colour. The distributions for Es and IfEs look very different 
and the KS test (Table [2]) quantifies that these distributions 
are not drawn from the same parent population. The isolated 
field elliptical galaxies show a bimodal distribution while Es 
show clearly only one peak. This result is also seen when 
the values of first quartiles of both distributions (Table [3| 
are compared. The value of the first quartile of our IfEs is 
almos t 0.3 magnitudes blu er than for the Es in agreement 
with lMarcum et al.l (|2004l ) who found evidence that 50 per 
cent of their sample of isolated early-type galaxies show blue 
global colours. 

The number density of isolated and non-isolated ellip- 
tical galaxies is quite similar for blue, B — R < 1.3, galaxies. 
The bimodality of the IfE distribution and the blue peak 
present in Fig. [6] is caused by the population of blue, faint 
and light IfEs. The colour distribution of IfEs suggests that 
some isolated field ellipticals have either formed later or have 
had merger activity at lower redshifts than non-isolated el- 
lipticals. We confirm whether this is the case in Section [S] 
where formation and merging times of IfEs are compared to 
non-isolated ellipticals. 

4.5 Dark Matter and Stellar Masses 

We continue our comparisons between IfEs and Es by study- 
ing their masses and composition. Both samples show very 
little cold gas: < 1.0 x 10 10 /i _1 Mq, with no ongoing star 
formation (the median star formation rate is ~ M©yr _1 ). 
This is typical for elliptical galaxies that are usually consid- 
ered as "red and dead" at z = 0. 

When dark matter halo masses are studied in detail, an 
interesting result is found. Isolated field elliptical galaxies 
reside mainly inside dark matter haloes that are lighter than 
7 x 10 12 /i _1 M s , while Es are found to reside more often 
in more massive haloes. The differences between the dark 
matter haloes of Es and IfE galaxies are striking, which is 
clearly visible in Fig. [7] where we plot the dark matter halo 
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Figure 7. Number density of dark matter halo mass for isolated 
field ellipticals (IfEs) and Es. 
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Figure 9. A scatter plot of stellar and dark matter halo mass. 
Yellow diamonds correspond to the sample of isolated field ellipti- 
cal galaxies (IfEs), while red circles correspond to Es. For clarity, 
only every 2nd galaxy has been plotted. 
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Figure 8. Number density of stellar mass for isolated field ellip- 
ticals (IfEs) and Es. 



mass distribution as a function of number density. It is no 
surprise that the KS test (Table [2} does not approve the null 
hypothesis, especially as the distribution of IfEs shows weak 
bimodality. The great difference in dark matter halo mass is 
interesting as the properties of galaxies are tightly related 
not only to the environment, but al so to the dark matter ha lo 
mass within the galaxy resides in (|Croton fc Farrarll2008h . 

While the dark matter haloes of IfEs are, in general, sig- 
nificantly lighter than haloes of Es, we see a different trend 
when stellar masses are compared. Even though the distri- 
butions in Fig. [8] may look somewhat similar, the KS test 
(Table [2]) disproves the null hypothesis with high probabil- 
ity. The differences in stellar mass distributions become very 
clear if we consider the quartile values of both samples. Iso- 
lated field elliptical galaxies have, in general, more stellar 
mass than Es (Table This is especially intriguing as IfEs 
reside inside rather light dark matter haloes (see Fig. [7|). 
To emphasize this difference we plot the dark matter halo 
mass versus the stellar mass for both samples in Fig. [9] IfEs 
are found to contain more stellar matter in respect to dark 
matter than Es. The stellar mass of IfEs grows almost lin- 



Table 2. Mean values of the properties and results of the 
Kolmogorov— Smirnov (KS) test: isolated field elliptical (IfEs) vs. 
control sample (Es) galaxies. 



Quantity 


IfEs 


Es 


KS Probability (p-value) 


Massj5M 


186.29 


6751.93 


< io- 10 


Massgy 


5.46 


3.92 


< IO" 10 


Mass CG 


0.32 


0.10 


< io- 10 


Age 


8.58 


9.32 


< io- 10 


Colour 


1.47 


1.58 


< io- 10 


M s 


-20.19 


-19.76 


< io- 10 


T 


-6.60 


-7.06 


0.39 



Note: Quantity Mass^i m refers to virial dark matter mass of the 

background halo, Massgr refers to stellar mass, while Massorj 
refers to the mass in cold gas. All masses are in units of 
10 10 /i — 1 Mq. Age is the mass- weighted age of a galaxy in units 
of 10 9 yr. Colour is B — R in absolute restframe (Vega) magni- 
tudes, and Ms is the absolute restframe (Vega) magnitude in the 
B— band (Buser B3 filter). T is the morphology, and it is the only 
quantity for which the null hypothesis of the KS test is approved. 



early with the dark matter mass with a mass (dark/stellar) 
ratio of ~ 4 x 10~ 2 . It is possible that the difference is re- 
lated to the formation and evolutionary paths of IfEs. This 
possibility is studied in detail in Section [5] 



4.6 Age distributions 

Distributions of mass weighted age (Fig. I10[) look rather 
similar, however, the KS test (Table [2} does not approve 
the null hypothesis. When quartile values (see Table [3} are 
studied it is obvious that IfEs have lower mass weighted age 
than Es, suggesting that IfEs are younger. Even IfEs seem 
to be slightly younger than Es in statistical sense, they cover 
roughly the same age range, only the extremes are missing. 
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Table 5. Results of the Kolmogorov— Smirnov (KS) test: isolated 
field ellipticals (IfEs) vs. the control sample (MEs) galaxies. 



Mass Weighted Age (10 9 yr) 



Figure 10. Number density of the mass weighted age for isolated 
field ellipticals (IfEs) and Es. 



Table 3. Quartile values of the properties of isolated (IfE) and 
control sample (Es) galaxies. 



Sample 


Quantity 


1 st Quartile 


Median 


3 rd Quartile 


IfEs 


Masss-M 


34.86 


120.58 


235.48 


Es 


Mass D M 


311.95 


1069.28 


4998.61 


IfEs 


Massgy 


1.69 


4.34 


7.94 


Es 


Massgy 


1.97 


2.84 


4.73 


IfEs 


MasscG 


0.10 


0.22 


0.44 


Es 


Massac 


0.04 


0.05 


0.09 


IfEs 


Age 


7.21 


8.84 


10.03 


Es 


Age 


8.44 


9.39 


10.41 


IfEs 


Colour 


1.28 


1.59 


1.64 


Es 


Colour 


1.56 


1.60 


1.63 


IfEs 


M s 


-20.70 


-20.18 


-19.64 


Es 


M s 


-20.14 


-19.65 


-19.28 


IfEs 


Ncomp 


3.00 


7.00 


15.00 



Note: MassoM is the virial dark matter mass of the background 

halo, Massgy - the stellar mass, Masspg - the mass in cold gas. 
All masses are in units of 10 10 h — 1 Mq. Age is the mass-weighted 
age of a galaxy in units of 10 9 yr. Colour is B — R in absolute 
restframe (Vega) magnitudes, and Ms is the absolute restframe 
(Vega) magnitude in the B— band (Buser B3 filter). Nqomp 
is the number of companion galaxies within the large (1.0/i — 1 
Mpc) avoidance sphere around the isolated field elliptical galaxy. 

5 FORMATION AND EVOLUTION OF IFES 

The analysis of simulated IfEs in the previous sections shows 
similarities in colour-magnitude diagrams and comparable 
ages, colours, and dark matter halo masses with observa- 
tional data. In general, we do not find big discrepancies 
between simulations and observations for the properties of 
IfEs we can compare. The KS tests between the IfEs and 
Es do not approve the null hypothesis, except for morphol- 
ogy. Therefore, IfEs selected with the criteria of Section [2. 2 1 
form a distinct class of objects and are significantly differ- 
ent from regular elliptical galaxies. As a next step in our 
analysis we use simulation data to study evolutionary paths 
and formation mechanisms of IfEs. Note that in this sec- 
tion, unless otherwise stated, we use the MEs sample as the 
control sample, which contains only those non-isolated ellip- 



Time 


D— value 


KS Probability (p— value) 


Za 


0.190 


5.7 x 10~ 8 


Zi 


0.239 


1.1 x 10" 5 


Z S 


0.353 


< io- 10 


Zl 


0.179 


6.5 x 10" 8 



Note: Quantities Xj , z a , zt and zi are the identity, assembly, for- 
mation and last merging times, respectively. 

tical galaxies that are the main galaxies of their dark matter 
haloes. 

5.1 Formation and Evolution Times 



iLucia fc Blaizotl (|2007l ) defined a set of times related to for- 
mation and evolution of dark matter haloes and galaxies 
that they reside in. We adopt the same definitions for con- 
venient and easy comparison, but we add one more quantity 
named last merging time, zi. Briefly, the different times are 
defined as follows: 

• Assembly time (z a ) is the redshift when 50 per cent of 
the final stellar mass is already present in a single galaxy of 
the merger tree. 

• Identity time (zi) is the redshift when the last major 
(the two progenitors both contain at least 20 per cent of the 
stellar mass of the descendant galaxy) merger occurred. 

• Formation time (z/) is the redshift when 50 per cent 
of the mass of the stars in the final galaxy at z — have 
already formed. 

• Last merging time (zi) is the redshift when the last 
merger occurred. 

We compute all four times, measured in redshifts, related to 
the formation and evolution of galaxies, and present numeri- 
cal results in Tabled Fig.[TT]shows cumulative distributions 
of the assembly (z a ), identity (zt), formation (z/), and last 
merging (zi) times. In Table [5] we show probabilities that 
the distributions of formation and evolution times of IfEs 
and MEs are drawn from the same parent distribution. The 
KS tests indicate great differences in all cases, implying that 
the formation and evolution of IfEs is different from MEs. 

Fig. [11] (top panel) shows that isolated field elliptical 
galaxies have assembled at lower redshifts than galaxies of 
MEs. Thus, in general, stars of IfEs form later than stars of 
MEs and thus IfEs are younger. This result is natural for 
hierarchical Cold Dark Matter models. If the conventional 
theory, that higher density areas collapse earlier, holds, this 
implies that IfEs have formed originally in less dense re- 
gions than MEs. Fig. 1111 (second panel) also shows that IfEs 
undergo their major merging events at significantly lower 
redshifts than the ellipticals of MEs. The difference in red- 
shift is significant (see Table [4} and implies that a different 
formation mechanism is behind the evolution of IfEs and 
MEs. 

Fig. [11] (third panel) shows that stars that will eventu- 
ally form an IfE galaxy are present already at higher red- 
shifts th an for MEs, in ag reement with observational find- 
ings (e.g. lReda et al.ll2005l ). These results indicate that IfEs 



12 



Sami-Matias Niemi et at 

Table 4. Statistics of formation and evolutionary times. 



Sample Time Mean Median 1 st Quartile 3 rd Quartile Mode Min Max Stdev 



IfEs 


Zi 


0.644 


0.408 


0.183 


0.687 


0.242 


0.064 


6.200 


0.836 


MEs 


Zi 


1.062 


0.624 


0.242 


1.386 


0.116 


0.020 


5.724 


1.089 


IfEs 


Za 


1.079 


0.989 


0.564 


1.504 


1.173 


0.020 


5.289 


0.710 


MEs 


Za 


1.303 


1.276 


0.755 


1.766 


1.276 


0.041 


4.520 


0.737 


IfEs 


Z f 


1.520 


1.386 


0.828 


1.912 


1.386 


0.012 


5.289 


0.831 


MEs 


Z f 


1.055 


0.687 


0.457 


1.386 


0.564 


0.041 


6.197 


0.893 


MEs* 


Z f 


3.214 


3.308 


2.831 


3.576 


3.576 


1.276 


6.197 


0.790 


IfEs 


Zl 


0.310 


0.208 


0.116 


0.408 


0.116 


0.020 


2.070 


0.325 


MEs 


Zl 


0.427 


0.230 


0.144 


0.509 


0.242 


0.020 


3.866 


0.483 



Note: Quantities Zi, z a , Zf and z; are the identity, assembly, formation and last merging times, respectively. 
MEs* refers to a sample where only haloes more massive than 5 X 10 12 /i — x Mq have been considered. 



Formation and Evolution Times 
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Figure 11. Cumulative distributions of assembly (z a ), identity 
(zi), formation (zf) and last merging (z;) times measured in red- 
shifts. The blue solid lines corresponds to IfEs, while the red 
dashed lines mark the MEs. 



can form stars more efficiently than MEs (see also Fig. I12[) . 
Note however, that we find the formation time (zf) to be 
highly de pendent on t he ma ss of the dark matter halo, as 
noted in iLucia et al l (|2006n . If we limit MEs to galaxies 
with dark matter halo mass greater than 5 x 10 12 /z _1 Mq we 
find the median formation time to be at very high redshift 
(z ~ 3.3). This shows that galaxies in massive dark matter 
haloes form the bulk of their stars already at very early cos- 
mic epochs, in agreement with the general hierarchical halo 
mass growth scenario. 



We define a major merger as an event where the two 
progenitors both contain at least 20 per cent of the stel- 
lar mass of the descendant galaxy. Almost half (~ 46.3 per 
cent) of isolated field elliptical galaxies have experienced a 
major merger at some point of their formation history. The 
percentage of major merging events for IfEs is higher than 
for MEs, where only one third (~ 33.3 per cent) of galaxies 
have experienced a major merging event. If all non-isolated 
E-type galaxies, independent of their luminosity or status 
in their dark matter haloes, are considered, only about four 
per cent experience a major merger. The difference is signif- 
icant and shows that it is possible to form elliptical galaxies, 
isolated or not, without a major merger via disc instabilities. 

Let us study the redshifts of the last merging event zi 
(sec Fig. [TT] bottom panel and Tables U and [SJ and the to- 
tal number of merging events. IfEs have their last merging 
event at lower redshifts than MEs. The median redshift of 
the last merging event is 0.21 and 0.28 for IfEs and MEs, 
respectively. This shows that IfEs have merging activity also 
at very late stages of t heir evolution in agree ment with 
iHernandez-Toledo et all l|200ct l. lHau et al.l (| 19991 ) quoted an 
age estimate of ~ 1.1 Gyr for isolated elliptical galaxies since 
the last merger. This age estimate from observational data 
seems to disagree with our median time of the last merging, 
which is more than twice as high: ~ 2.5 Gyr. However, we 
do find IfEs that have experienced their last merging event 
only ~ 0.3 Gyr ago. Thus, the discr epancy can arise from 
the small number of observed IfEs in lHau et al.l l|l999l ). 

It is possible that evolution of IfEs is suppressed by the 
low density environment they reside in, thus the major merg- 
ing events happen later stage in their evolution if at all. The 
larger fraction of IfEs having major mergers than MEs can 
also be due to the requirement of the magnitude gap not only 
the region they formed in. IfEs of denser areas with compa- 
rable sized galaxies must clean their environment effectively, 
leading to a larger fraction of major merging events. 

While studying the redshifts of the last merging events 
we noticed that 6 IfEs did not experience a single merg- 
ing event during their evolution. These galaxies have de- 
veloped in truly isolated areas; although, maybe, a greater 
mass resolution in simulations would reveal one or more mi- 
nor merging events. Despite the resolution limitations, the 
result is in agreement with observations where some IfEs 
do not show any signs of merging activity. Typically, simu- 
lated IfEs experience one to 10 merging events (above the 
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Figure 12. Mass assembly of isolated field ellipticals (IfEs) and 
MEs as a function of redshift z. Lines show the median value at 
given redshift. 

mass limit of 1O s /i _1 M0) during their lifetime, the average 
being 7.6 and the median six mergers. These numbers are 
significantly higher than for MEs, for which we find the av- 
erage and median of 4.6 and one merging event, respectively. 
This suggests that IfEs accrete multiple smaller haloes (and 
galaxies) during their formation, cleaning their neighbour- 
hood from dwarf galaxies effectively. This leads to a question 
how effectively IfEs accrete mass as a function of time com- 
pared to MEs. 

5.2 Mass Assembly 

Fig. [12] shows the dark and stellar matter mass assembly 
as a function of redshift for both IfEs and MEs. The figure 
shows that IfEs start to form around the same epoch as 
MEs, however, IfEs accumulate stellar matter much faster. 
According to our findings IfEs can form stars more efficiently 
than MEs, while MEs seem to accrete dark matter slightly 
faster than IfEs. We confirm th at IfEs form t he bu lk of their 
stars at z > 2, as suggested in iReda et al.l (|2005r ). We also 
note that at redshift one IfEs have formed over half of their 
stars (stellar mass) and have gathered as much as 80 per 
cent of their final dark matter. The galaxies of MEs have 
accreted roughly the same fraction of dark matter as IfEs 
at z ~ 1, however, they have formed as little as 30 per cent 
of their stars compared to the final stellar matter at z = 0. 
This difference is significant and shows that IfEs form their 
stars quickly and are rather dark-matter poor compared to 
other elliptical galaxies (see also Fig. [5] for stellar vs. dark 
matter). It is also noteworthy that IfEs continue to accrete 
dark matter till z — while MEs have gathered ~ 99 per 
cent of their final dark matter already at z ~ 0.5. All this 
points towards a different formation mechanism for these 
two galaxy classes and suggests that late merging events 
can be a significant part of the evolution of an IfE. 

To further illustrate the mass assembly of IfEs, we 
plot the sample of blue, light and faint IfEs separately in 
Fig. 1131 We have divided the IfEs sample into two: light 
(M DM < lO 12 /i _1 M ) and heavy {M DM > 1O 12 /i _1 M ) 
IfEs. Fig [13] shows that there is a significant difference in 
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Figure 13. Mass assembly of isolated field ellipticals (IfEs) as a 
function of redshift z when the sample of IfEs has been divided 
by dark matter halo mass. Lines show the median value at a given 
redshift. 



mass accretion of light and heavy IfEs. The stellar mass of 
heavy IfEs follows closely the dark matter mass accretion, 
while light IfEs evolve differently. The heavy IfEs have cre- 
ated half of their stellar matter already at z ~ 1.6, while 
the light IfEs have created barely 10 per cent of their stella r 
matter, in agreement with the findings of lTreu et al.l {2005). 
We note that the light IfEs create half of their stellar matter 
by z ~ 0.7. The heavy IfEs also form stars extremely effi- 
ciently compared to the light IfEs, as the stellar mass follows 
dark matter accretion closely. Fig. [T3] also shows that heavy 
IfEs host older stellar populations than light IfEs, as ~ 99 
per cent of their stellar mass has been accumulated already 
by z ~ 0.3. Thus, more massive (and luminous) IfEs have 
old stellar populations, while the lighter ones have formed a 
signification fraction of their stellar mass relatively recently. 



5.3 Merger Histories of IfEs 

So far we have shown that minor and major merging events 
take place during the formation of an IfE. Thus, in the 
following subsections, we try to identify general formation 
mechanisms based on different types of merger trees and 
merging events. We study merger trees of all IfEs individu- 
ally and identify three typical formation scenarios. 

We illustrate merger histories in merger plots. The iso- 
lated field elliptical itself lies at the bottom of the plot at z = 
0, and all its progenitors (more massive than 10 8 /i _1 Mq) are 
plotted upwards going back in time. Galaxies with stellar 
mass larger than 10 9 /i _1 Mq are shown as symbols, and are 
colour-coded as a function of their rest-frame B — R colour. 



5.3.1 Solitude 

Fig. ll4l shows an example of an isolated field elliptical galaxy 
that has developed undisturbed, alone, and has not under- 
gone even a single merging event. We group IfEs that form in 
this fashion to a single group that we call solitude formation 
class. 
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Even though the Fig. [14] does not show a single merg- 
ing event during the formation history, this may not be the 
whole story; the plot does not show mergers smaller than the 
resolution limit (~ 10 8 /i _1 Mq) of the Millennium Simula- 
tion. Therefore, it is possible or even likely that IfEs belong- 
ing to this class have had one or even several minor merging 
events during their evolution. Even so, the merging events 
would have involved very light dark matter clumps, and it is 
likely that no significant observational evidence would exist. 
In this context small merging events can be interpreted as 
accretion of dark matter making solitude a proper formation 
mechanism. 

We identify six IfEs, corresponding to mere two per cent 
of all IfEs, that belong to the solitude class. This implies that 
either IfEs that develop truly alone in underdense regions are 
extremely rare or that some of the IfEs have been misclas- 
sified. It is possible that some IfEs with one or more small 
mergers could belong to this formation class, especially if 
the motivation of classification is based on observability of 
each class. It is not obvious how massive a merging event is 
required to find observational evidence of a merger compli- 
cating the classification. 

The six IfEs identified all show unsteady evolution in 
their colour as a function of redshift. This is somewhat 
surprising as one would assume that a passively evolving 
galaxy should show a steady colour evolution from blue to 
red while the stellar population ages. The bulge formation 
of solitude class is assumed to happen via disc instabilities. 
Solitude IfEs reside inside lighter dark matter haloes than 
IfEs of other class, with typical dark matter halo mass of 
or even less. This makes IfEs that be- 
long to the solitude class the lightest IfEs in the MS. The 
formation, when the dark matter mass grows larger than 
1O 8 /i -1 M0, epoch of solitary IfEs is in the redshift range 
3.6 < 2 < 5.3. However, we find one solitude IfE that formed 
as late as z ~ 2.5, making it a very late bloomer, show- 
ing that, IfEs can have formed very rece ntly. Solitude IfE s 
are in agreement with observations of e.g. lAars et all (|200lT ) 
and lDenicolo et al l (|2005h . who did not find any evidence of 
merging activity while studying their samples of IfEs. 

5.3.2 Coupling 

Fig. [15] shows an example of an isolated field elliptical that 
has undergone at least one "equal" size merger during its 
evolution. IfEs that have experienced at least one equal sized 
merger comprise a formation class named coupling. 

Our definition for an equal size merging is the following: 
the lighter of the merged galaxies had at least 50 per cent of 
the stellar mass of the more massive one. This is far larger 
than in our definition of a major merger (20 per cent) that 
was applied when calculating the identity time Zi. Thus, an 
equal sized merger guarantees that the merging event has 
had a great impact not only on the morphology, but also on 
the evolution and other properties of the IfE. 

The motivation for this formation class resides in ob- 
servations; a major merging event should be visible in ob- 
servational data. Thus, an equal sized merging should leave 
distinct marks to the descendant galaxy, which should be 
observable for at least few Gyrs if not longer. It has been 
suggested that isolated field ellipticals have formed in rel- 
atively recent mergers of spiral galaxy pairs, i.e. in merg- 
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Figure 14. Example of a merger tree of an isolated field elliptical 
galaxy that has developed alone without any significant merging 
events. Symbols are colour-coded as a function of the B — R colour 
and their area scales with the stellar mass. Only progenitors more 
massive than 10 8 /i _1 Mn are shown. 



ing of two c o mpar able si zed galaxies. | Marcum et alJ l]2004h . 
iReda et all (|2007l l and iKautsch et al l (|200Sf ) have found 
several isolated galaxies that show signs of re cent mor- 
phological disturban ces, while lHau et ahl (|l999h (see also 
lHau fc Forbes! 120061 ) have found that ~ 40 per cent of 
isolated galaxies show kinematically distinct cores (KDC). 
KDCs are generally be lieved to be the result of a major or 
an equal sized merger (jrlernquist fe Barneslll99ll 'l. Thus, it 
is possible that these observations have already identified 
galaxies that belong to this formation class. However, it is 
also possible that KDCs can form in an early colla pse with- 
out subsequent mergers (|Harsoula Sz Voglislll99ct ). compli- 
cating the identification of the formation mechanism. 

We identify in total 93 IfEs, corresponding to ~ 32 per 
cent of all IfEs, that belong to the coupling class of formation 
scenarios. IfEs that belong to the coupling formation class 
show colour evolution that is in agreement with conventional 
theory. The time from the last equal sized merging correlates 
well with the colour of the galaxy at z = 0; galaxies with 
late merging are bluer than galaxies that experienced their 
equal sized merging a long time ago. The mean redshift of 
the last equal sized merger is 1.09 while the median is 0.62. 
IfEs that belong to this formation class start to form usually 
in the redshift range 3.6 < z < 5.3. We also identify a few 
IfEs that have started to form as early as z ~ 8, making 
IfEs that belong to the coupling class older, in a statistical 
sense, than IfEs that belong to the solitude formation class. 

The coupling formation mechanism can explain sev- 
eral observed IfEs. The colour-magnitude diagrams of 
IReda et al.l (|2004f ) show a slope and scatter that is in agree- 
ment with equal-mass mergers. They find 11 per cent of their 
IfEs to show bo xy isophotes th a t can form when equal sized 
galaxies merge. M arcum et al l (|2004l ) argue that all except 
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Figure 15. Example of a merger tree of an isolated field elliptical 
galaxy that has undergone an equal sized merger. Symbols are 
colour-coded as a function of the B — R colour and their area 
scales with the stellar mass. Only progenitors more massive than 
10 8 /i _1 M Q are shown. 
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Figure 16. Example of a merger tree of an isolated field elliptical 
galaxy that has undergone multiple merging events, but not any 
equal sized ones. Symbols are colour-coded as a function of the 
B — R, colour and their area scales with the stellar mass. Only 
progenitors more massive than 10 8 /i -1 Mq are shown. 



one of their IfE have luminosities that would, at most, be 
consistent with a single equal- mass merger event. Note that 
all these observations are best explained by the IfEs of the 
coupling class. 

5.3.3 Cannibalism 

Fig. [16] shows a typical isolated field elliptical galaxy that 
has developed and accreted dark matter and stellar mass 
through multiple, small and larger merging events, but has 
not experienced any equal sized mergers. All IfEs that form 
via multiple mergers form a formation class called cannibal- 
ism. 

The merger trees of cannibal IfEs show a large number 
of mergers, with also a relatively large ones. It is obvious 
from the merger trees that the colour of the IfE does not 
change due to a minor merger. This is due to the fact that 
many small merging events are likely to be dry, and there- 
fore, do not induce significant star formation that would 
make the global colour of the IfE bluer. However, the largest 
merging events can have a significant impact on the IfE's 
morphology and colour (see Fig. [TS] and the merging event 
around 0.8 Gyr ago). Some merging events, visible in the 
merger trees of cannibal IfEs, can be classified as major 
mergers, and these could be visible in observational data. 

Most (194) of the IfEs in the MS belong to the cannibal 
class, corresponding to ~ 66 per cent of all IfEs. The large 
number of IfEs belonging to this class shows that merging 
events are important for the formation and evolution of iso- 
lated field elliptical galaxies. Cannibals often show a large 
number of merging events, more than 20, while the number 
of major merging events ranges from one to three. The pos- 
sibility of major mergers can complicate the identification of 



IfEs of this class, especially as we cannot identify any pre- 
ferred time for the last major merger. In general, IfEs of the 
cannibal class form early, in the redshift range 5 < z < 12. 
However, we identify few cannibal IfEs that form as late as 
2 ~ 2.4. The cannibal IfEs that form late could also be clas- 
sified as solitudes, especially as these IfEs do not, in general, 
show a single major merging event only few minor ones. 

This formatio n class can explain a few observed IfEs. 
iReda et all \2004 ) found low-luminosity dwarfs close to an 
isolated galaxy that avoid accretion. Three of their IfEs show 
high central space density, but they also found quite strong 
morphological disturbances requiring an accretion of a fairly 
large galaxy. This is in agreement with cannibal IfEs that 
have had a major merging event that would explain the 
morphological disturbances, while a large number of smaller 
satellite dwarfs could have survived due to the large dark 
matter halo and long dynamical friction times. 



6 DISCUSSION 

6.1 Population of Blue and Faint IfEs 

Figs. [T] and [2] show a separate population of blue and faint 
IfEs. This population was also confirmed to comprise iso- 
lated field elliptical galaxies that reside in light, Mum < 
1 12 ft~ 1 Mp), dark m atter haloes. Moreover, th e IfE samples 
of lReda et all (|2004T ) and lMarcum et all (|2004f) do not reveal 
a single observed IfE that would belong to this predicited 

pop ulation of IfEs. 

iHernandez- Toledo et all |2008) studied isolated galax- 
ies and used modern observations of Sloan Digital Sky Sur- 
vey (SDSS) Data Release (DR) 6. The base sample of their 
galaxies were the Catalogue of Isolated Galaxies (CIG) 
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in the Northern Hemisphere combiled by iKarachentseval 
J 19731). Although th e isolation criteria in the study of 
IKarachentseval l|l973l ) and in our study are not the same, 
it is still interesting to compare if any of the elliptical galax- 
ies could belong to t he population of faint and blue If Es that 
our results predict. IHernandez- Toledo et al.l (|2008T ) noted 
that four early-type galaxies of the CIG sample showed blue 
(g — i< 1.0) colours. Three of their blue galaxies were found 
to be disky while the other one was found to be boxy. From 
the morphological inspections they concluded that three of 
the blue galaxies were ellipticals while one was SO-type. 
Thus, below we will inspect if any of the three bl ue elliptical 
galaxies noted in IHernandez- Toledo et all f2008), could be- 
long to the faint and blue population of IfEs that simulations 
predict. 

As the magnitudes of IHernandez- Toledo et all (2008) 
are in the SDSS system, we need a method to transform their 
magnitudes for a compar ison. For a crude com parison we 
can use galaxy colours of iFukugita et al.l (|l995l ) and derive 
a transformations: 



6.2 Observing a Formation Class 



Me 



M r - 0.35 



(2) 



for R— band magnitudes for a typical elliptical. Because of 
our IfEs show bluer colours than typical ellipticals, this 
transformation is not very accurate, however, for our pur- 
poses it should be adequate. To transform the SDSS colours, 
we use the following equation: 



B-R, 



-i) + 0.44 



(3) 



which has been derived from the works of IFukugita et al.l 
( 1996) an d l Jester et all (|2005l ). Again, we note that this may 
not provide exact colour transformation, but should provide 
satisfactory conversion for our crude comparisons. 

The three isolated ellipticals of lHernandez- Toledo et all 

(2008) show g - i colours of 0.91, 0.98 and 1.02. With the 
equation [3] we can approximate that these colours corre- 
spond to B — R colours of 1.35, 1.42, and 1.46. Clearly at 
least one of their galaxy is blue enough (B — R ^ 1.4) to be- 
long to the separate population of faint and blue IfEs. Given 
the inaccuracy of equation [3] it is possible that all three of 
their galaxies would belong to the blue population o f IfEs. 
The blue galaxies of IHernandez- Toledo et al.l ((2008) have 
absolute r— band magnitudes —20.84, —21.24, and —16.73, 
respectively. With the help of equation [2] we can approxi- 
mate their R— band magnitudes to be —21.19, —21.59, and 
— 17.08, respectively. Thus, at least two of their galaxies are 
faint enough {Mr > —21.5) to be part of the population. 

Above magnitudes and colours show that possibly at 
least one and up to three isolated field elliptical galaxies that 
belong to the predicited population of faint and blue IfEs 
has already been observ e d. On e of the isolated ellipticals of 
IHernandez- Toledo et al.l (|2008l ) is significantly fainter than 
any of our IfEs, and would not qualify as an IfE in our 
study. Moreover, as the isolation criteria of CIG galaxies 
are different that ours, we cannot conclude without a doubt 
that any of the blue and faint population IfEs have been 
observed. Thus, more observations of blue and faint isolated 
elliptical galaxies are required to confirm our prediction. 



Section \5. 31 introduced three typical formation classes of iso- 
lated field elliptical galaxies; solitude, coupling, and canni- 
balism. Our results show that the majority of IfEs experience 
numerous merging events during their evolution and belong 
to the cannibalism class. A smaller, yet a significant frac- 
tion of IfEs belong to the coupling class, whose members 
have experienced an equal sized merging, while only a small 
fraction of IfEs show insignificant merging activity and be- 
long to the solitude class. If merging events are so important 
for the evolution of IfEs, an obvious question remains: Is it 
possible to identify the formation scenario based on obser- 
vational data? Therefore, in this section we briefly discuss 
the properties of IfEs that can be observed and that at the 
same time would readily indicate the formation mechanism 
of the IfE. 

IfEs that belong to the coupling class form later than 
cannibal IfEs and are found to reside in lighter dark matter 
haloes than cannibal IfEs: the median masses are 54 and 
200 x lO lo /i _1 M0, respectively. The coupling class IfEs can 
therefore be merger remnants of nearby galaxy pairs having 
a modest sized dark matter halo at redshift zero. Unfortu- 
nately, dark matter haloes cannot be observed directly. IfEs 
of equal sized mergings can show a weak X-ray emission; 
however, it is not clear how well this could be detected if 
at all. Moreover, differentiating coupled IfEs from cannibals 
might be difficult from X-ray data only. 

The mean number of companion galaxies of IfEs that 
have formed through coupling is lower than for the cannibal 
IfEs, but larger than for the solitary IfEs. Unfortunately, the 
distribution of stellar mass, B— band magnitude and colour 
of IfEs of the coupling class is indistinguishable from the 
cannibal IfEs. The redshifts of equal sized merging events 
suggest that to be able to identify galaxies belonging to the 
coupling class, observations should concentrate on interme- 
diate redshifts in the range 0.25 < z < 1.4. However, to find 
evidence of equal sized merging may not be simple; it is not 
clear how different are the traces an equal size merger leaves, 
compared to a major merger. Study of stellar populations 
could provide a way to identify IfEs that have experienced 
an equal sized merging, but this may not be applicable for 
high redshifts. 

A cannibal IfE can be a merger remnant of a small 
or compact group. We find cannibal IfEs to reside in 
more massive and larger dark matter haloes than other 
types of IfEs. However, their dark matter haloes are 
less massive than haloes of typical groups, leaving only 
small or compact groups to consider. Multiple merging 
events of cannibal IfEs could also be visible in their 
stellar populations, especially if merging events were gas 
rich, so-called wet mergers. Observational evidence shows 
that it i s possible to detect IfEs in X-ray observa- 
tions (e.g. iMulchaev fe Zabludoflj|l999l; ISivakoff et al.ll2004l : 



lO'Sullivan fe Ponmanl 120041 ; iMemola et alj 120091 ). Due to 

their more massive dark matter haloes, X-ray bright IfEs 
are the best candidates for the cannibal formation class. 



6.3 Comparison to Fossil Groups 

Have fossil groups formed in a similar way as isolated field 
elliptical galaxies? Are IfEs an intermediate product while 
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they develop to become fossil groups, or vice versa? These 
questions are justified as the selection criteria of IfEs are 
very similar to those used for the identification of fossil 
groups, especially if only optical data are available. Objects 
of both classes are required to show a large magnitude gap 
between the brightest and the second brightest galaxy in 
optical. In general, fossil groups are also required to show 
extended X-ray emission while it is not demanded for IfEs. 
However, several IfEs have shown some amount of extended 
X-ray emission even it is not required. Due to similar se- 
lection criteria objects of both classes might have similar 
evolution and assembly histories and a similar formation 
mechanism, being actually the same class seen only at dif- 
ferent phases of evolution. 

The majority of fossil groups seem to have experi- 
enced the last major merging event longer than 3 Gyr 
ago and th ey have assembled half of their final mass 
by z ^ 0.8 dvon Benda-Beckmann et aDl2008l ). Moreover, 
Ivon Benda-Beckmann et al. |2008l ) found from simulations 
that only 15 per cent of fossil groups experience the last 
major merger less than 2 Gyr ago, and at least 50 per cent 
had the last major merger longer than 6 Gyr ago. These 
timescales are in modest agreement with our findings. Our 
results show that IfEs have experienced their last major 
merging, on average, ~ 6 Gyr ago, while IfEs have half of 
their final mass assembled by z ~ 1. For f ossil groups this 
num ber is ~ 0.6 ( Di'az-Gimc nez et alJ l2008) . 

iDfaz-Gimenez et al.l (|2008l ) study the evolution of the 
brightest galaxies of fossil groups. They identified fossil 
groups from the Millen nium Simulation and a dopted the 
same galaxy catalogue l|Lucia fc Blaizotl l2007h as in this 
study, making it interesting to compare their findings to 
ours. When comparing the medians of assembly, formation 
and identity times we find that isolated field elliptical galax- 
ies have assembled their stars earlier than the brightest 
galaxies of fossil groups; z a ~ 1.1 and ~ 0.6, respectively. 
However, fossil groups form significantly earlier than IfEs; 
Zf ~ 3.6 and ~ 1.5, respectively. The median time of the 
last major merging of the IfEs is almost twice (zi ~ 0.6) 
the identity time of the brightest galaxies of fossil groups 
(zi ~ 0.3). These discrepancies in formation and evolution- 
ary times cast a serious doubt over the idea of similar for- 
mat ion mechanisms. 

Ivon Benda-Beckmann et alj |2008l ) argue that the pri- 
mary driver for the large magnitude gap is the early infall of 
massive satellites that is related to the early formation time 
of fossil groups. We do not find infall of massive satellites at 
early time in our sample of simulated isolated field elliptical 
galaxies. The early evolution of an IfE includes relatively mi- 
nor mergers, except for some IfEs that experience an equal 
sized merging. However, the time of the equal sized mergers 
is usually in a later stage of the development of the IfE, not 
early as argued for fossil groups. Our study shows that the 
primary driver for the large magnitude gap of IfEs is either a 
merging of a comparable sized galaxy pair or effective mass 
accretion and sweeping up of surrounding galaxies (coupling 
and cannibalism, respectively). 

iDariush et all (120071 ) argue that fossil groups can be 
identified from dark matter only simulations if one selects 
dark matter haloes more massive than 5x W 13 Ii~ 1 Mq. They 
argue that above that limit, all optical fossil groups in the 
Millennium Simulation have enough hot gas to show X-ray 



emission, therefore qualifying as X-ray fossils. Our Table [3] 
shows that the median mass of the dark matter haloes of 
IfEs is only ~ 1.2 x lO 12 ft _1 M , while even the 3rd quar- 
tile is just ~ 2.4 x W 12 h~ 1 MQ. This comparison shows that 
IfEs reside in significantly lighter dar k matter haloes than 
the brightest galaxies of fossil groups. IDariush et al.l (|2007l ) 
argued that fossil groups have formed early and that more 
than ~ 80 per cent of their mass accumulated as early as 
4 Gyr ago. Moreover, they suggest that X-ray fossil groups 
are not a distinct class of objects but rather that they are 
extreme examples of groups which collapse early and expe- 
rience little recent growth. This formation scenario does not 
apply for IfEs, as cannibal IfEs can have merging activity 
till the redshift z ~ 0, and IfEs start to form later than fossil 
groups. 

Although a large magnitude gap seen in both fossils 
and IfEs could imply a similar formation mechanism, the 
above comparison does not support this idea. Considering 
the assembly and formation times it seems unlikely that 
fossil groups, and especially their brightest galaxies, could 
share the same formation mechanism as isolated field el- 
liptical galaxies. The large differences in dark matter halo 
masses of fossil groups and the most massive IfEs makes it 
improbable that fossil groups and IfEs are the same class of 
objects seen at different phases of their evolution. We there- 
fore conclude that fossil groups and IfEs form two distinct 
classes. 



7 SUMMARY AND CONCLUSIONS 

The aim of this paper is twofold: 1) to compare simulated 
field elliptical galaxies with observed ones and to make pre- 
dictions on their properties, and 2) to define the formation 
mechanism and history of isolated field elliptical galaxies. 
We also discuss if isolated field elliptical galaxies are related 
to fossil groups and if they can share a common formation 
mechanism. 

Our results show that it is possible to identify isolated 
field elliptical galaxies (IfEs) from cosmological TV— body 
simulations with semi-analytical models of galaxy forma- 
tion, that have similar properties to observed IfEs. We show 
that simulated IfEs are in good agreement with observa- 
tions when similar identification criteria are adopted. The 
colour-magnitude diagram of simulated IfEs agree with ob- 
servations, and the average colour of simulated and observed 
IfEs are within their error limits when all simulated IfEs are 
considered. Unfortunately, observational data sets are still 
small complicating more detailed comparison. An agreement 
in age and mass estimates of IfEs between observations and 
simulations is found. However, the age comparisons are less 
robust due to different age definitions. 

Our results show that isolated field ellipticals are very 
rare objects; we find their total number density to be as 
low as ~ 8.0 x 10 _6 h 3 Mpc -3 . Our result agrees with ob- 
servational estimates of number of IfEs, which however, are 
inaccurate at best. Our IfEs have a small number of com- 
panion galaxies, ranging from only a few dwarf companions 
to as much as about 20. Thus, IfEs are not completely iso- 
lated, although they are found to be located in underdense 
regions. Our results show that IfEs reside in relatively light 
dark matter haloes. However, at the same time, the baryonic 
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to dark matter ratio is higher in IfEs than in Es generally. 
The stellar mass of IfEs grows almost linearly with the dark 
matter mass with a mass ratio (dark/stellar) of ~ 4 x 10 -2 , 
whereas our comparison ellipticals have a lower stellar to 
dark matter ratio. Therefore, IfEs are good candidates for 
galaxies with low dark matter mass to stellar light ratio. We 
also find a flat B— band luminosity function for our simu- 
lated IfEs. 

When studying the basic properties of IfEs we find 
that IfEs populate different regions in colour-magnitude 
diagrams than regular elliptical galaxies, which are found 
mostly in the red sequence. This is due to the bimodality of 
the distribution of colours, as IfEs populate not only the red 
sequence, but also the blue cloud. From simulation data we 
find a separate population of blue and faint IfEs. On average 
IfEs are found to be bluer than our control sample ellipticals, 
howevever, this result is biased because of the separate pop- 
ulation of IfEs. The bluest IfEs are found to reside in light 
dark matter haloes, while red IfEs are usually found inside 
more massive (Mom Z 10 12 hT 1 Mq) dark matter haloes. 
We note that simulations predict a previously unobserved 
population of blue, dim and light galaxies that fulfill ob- 
servational criteria to be classified as isolated field elliptical 
galaxies. These galaxies have only a few companions, which 
are usually located many times further away than the virial 
radius of their dark matter halo. These blue, dim and light 
IfEs have formed their stars only lately and have continued 
to accrete dark matter mass till redshift zero. This popula- 
tion of IfEs is very interesting as it has not been detected 
yet in observations. 

Our results also show that IfEs start to form around 
the same epoch as the galaxies of the control sample (MEs); 
however, IfEs seem to accumulate stellar mass much faster. 
IfEs form stars more efficiently than MEs, while MEs ac- 
crete dark matter slightly faster than IfEs. We confirm that 
IfEs form the bu lk of their stars at z > 2, as suggested in 
iRedaet al.ll|2005h . By the redshift one IfEs have formed over 
half of their stars (stellar mass) and have gathered as much 
as 80 per cent of their final dark matter. The galaxies of MEs 
sample have accreted roughly the same fraction of dark mat- 
ter as IfEs by z ~ 1. However, they have formed as little as 
30 per cent of their stars compared to the final stellar mat- 
ter at z — 0. This difference shows that IfEs form their stars 
quickly and are rather dark-matter poor compared to other 
elliptical galaxies. Moreover, IfEs continue to accrete dark 
matter till z — while MEs have gathered ~ 99 per cent 
of their final dark matter already by z ~ 0.5. We note that 
more massive (and luminous) IfEs have older stellar pop- 
ulations, while the lighter ones have formed a signification 
fraction of their stellar mass relatively recently. 

While studying the evolution of IfEs we note that al- 
most half (~ 46 per cent) of IfEs have experienced at least 
one major merger during their formation history, while only 
about four per cent (Es) up to one third (MEs) of control 
sample ellipticals experience a major merger. Major merging 
events happen later in IfEs evolution than for control sam- 
ple galaxies, the average of the latest major merging being 
z ~ 0.6 for IfEs while it is z ~ 1.1 for MEs. We also find 
IfEs that have not experienced a single merger event above 
the mass resolution limit of the MS during their evolution. 
Therefore, it is possible to form elliptical galaxies without 
major mergers. 



When inspecting the merger trees of simulated IfEs, we 
identify three typical formation scenarios: solitude, coupling, 
and cannibalism, which can all lead to a formation of an IfE. 
The scenarios range from a solitary growth (solitude class) 
with quiet mass accretion and star formation to more violent 
evolution with multiple mergers (cannibalism class). We also 
identify a formation scenario where two comparable sized 
galaxies merge to form an IfE (coupling class). Our merger 
trees show that merging events are important for IfEs that 
form through cannibalism or coupling. All three formation 
classes are in agreement with observational findings. 

Comparison between isolated field elliptical galaxies 
and fossil groups show that these two classes are distinct. 
Galaxies of both classes show some similarities, but many 
properties and evolution times are significantly different. 
IfEs reside in significantly lighter dark matter haloes and we 
do not find an infall of massiv e galaxies at early times for 
IfEs as argued for fossil groups l|von Benda-Beckmann et al.l 
2008). However, we cannot exclude that some fossil groups 
could not share the formation mechanism of the most mas- 
sive IfEs, namely cannibalism. 
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